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This invention was made with Government support under NIH Grant 
Nos. DK-47324, HL-47401, HL-67647 and HL-68658. The Government has 
certain rights in the invention. 

BACKGROUND OF THE INVENTION 

10 1 . Field of the Invention 

This invention relates to methods and systems for measuring 
mechanical property of a vascular wall and a method and system for determining 
health of a vascular structure. 

2. Background Art 

15 The following references are referenced herein: 

A.J. Bank et al., "In Vivo Human Brachial Artery Elastic Mechanics 
Effects of Smooth Muscle Relaxation," 1999; Circulation 100:41-47; 

D.H. Bergel, "The Static Elastic Properties of the Arterial Wall," J. 
Physiol., 1961; 156:445-457; 
20 C. Bilato et al., "Atherosclerosis and Vascular Biology of Aging," 

Aging (Milano) 1996; 8:221-234; 

O. Bonnefous et al., "Non Invasive Echographic Techniques for 
Arterial Wall Characterization," IEEE Ultrasonic Symposium, .1996: 1059-1064; 

J. Blacher et al., "Carotidarterial Stiffness as a Predictor of 
25 Cardiovascular and All-Cause Mortality in End-Stage Renal Disease," 
Hypertension, 1998; 32:570-574; 
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J. Blacher et al., "Impact of Aortic Stiffness on Survival in End-Stage 
Renal Disease," Circulation 1999; 99:2434-2439; 

O. Bonnefous et al., "New TDI Developments for Vascular and 
Cardiac Applications," IEEE Ultrasonic Symposium, 2000: 1285-1290; 
5 A. Bruel et al., "Changes in Biomechanical Properties, Composition 

of Collagen and Elastin, and Advanced Glycation Endproducts of the Rat Aorta in 
Relation to Age," Atherosclerosis 1996; 127:155-165; 

D. Buprez et al., "Relationship Between Periventricular or Deep 
White Matter Lesions and Arterial Elasticity Indices in Very Old People," Age and 
10 Ageing, 2001; 30:325-330; 

A. Eriksson et al., "Arterial Pulse Wave Velocity with Tissue 
Doppler Imaging," Ultrasound in Med. & Biol., 2002; Vol. 28, No. 5:571-580; 

G. Faury, "Function-Structure Relationship of Elastic Arteries in 
Evolution: From Microfibrils to Elastin and Elastic Fibres," Pathol. Biol., 2001; 
15 49:310-325; 

Y.C. Fung, "Biomechanics: Mechanical Properties of Living 
Tissues," 2nd Ed., Spring- Verlag, 1993: 321-391; 

G. Guerin et al., "Arterial Stiffening and Vascular Calcifications in 
End-Stage Renal Disease," Nephro Dial Transplantation, 2000; 15:1014-1021; 
20 Hardung V., "Propagation of Pulse Waves in Visco-Elastic Tubing," 

American Physiological Society, Handbook of Physiology, Section 2, 
Circulation, 1962, Vol. 1, eds., Hamilton, W.F. and Dow, P., 107; 

D.R. Kaiser et al., "Brachial Artery Elastic Mechanics in Patients 
with Heart Failure," 2001; Hypertension 38:1440-1445; 
25 K. Konner et al., "The Arteriovenous Fistula," J. Am. Soc. 

Nephrol., 2003; 14(6): 1669-80; 

G.J. Langewouters et al., "The Static Elastic Properties of 45 Human 
Thoracic and 20 Abdominal Aortas In Vitro and the Parameters of a New Model," 
J. BlOMECH., 1984; 17-425-435; 
30 M.A. Lubinski et al., "Speckle Tracking Methods for Ultrasonic 

Elasticity Imaging Using Short Time Correlation," IEEE Trans. Ultrason., 
Ferroelect., Freq. Contr., 1999, Vol. 46, pp. 82-96; 
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A.J. Luik et al., "Arterial Compliance in Patients on Long- 
Treatment-Time Dialysis," Nephrol. Dial Transplant, 1997; 12:2629-2632; 

J.J. Mai et al., "Strain Imaging of Internal Deformation," 
Ultrasound in Med. & Biol., 2002; Vol. 28, Nos. 11/12:1475-1484; 
5 M. Persson et al. , "Estimation of Arterial Pulse Wave Velocity With 

A New Improved Tissue Doppler Method," Proceeding of the 23rd Annual 
EMBS International Conference, 2001:188-191; 

H. Taniwaki et al., "Femoral Artery Wall Thickness and Stiffness in 
Evaluation of Peripheral Vascular Disease in Type 2 Diabetes Mellitus," 
10 Atherosclerosis, 2001; 158:207-214; and 

S. Timoshenko et al., "Theory of Elasticity," 3rd Ed., McGraw 
Hill, New York, 1970. 

Arterial compliance has been shown to be a strong indicator of 
vascular disease; cardiovascular disease, peripheral vascular occlusive disease, 

15 diabetes, and renal failure. Changes in the ratio of collagen to elastin in the 
extracellular matrix of the arterial media is believed to be one of the causes of 
arterial stiffness (Faury 2001; Bilato and Crow 1996; Bruel and Oxlund 1996). By 
measuring mechanical properties of tissue, elasticity imaging could non-invasively 
monitor vascular pathologies developing within the vascular wall. Previous attempts 

20 at non-invasive vascular elastic imaging include arterial wall motion estimation 
(Bonnefous et al., 1996; Taniwaki et al., 2001; Luik et al., 1997; Guerin et al., 
2000), intraparietal strain imaging (Bonnefous et al., 2000) and pulse wave velocity 
measurement (Eriksson et al., 2002; Persson et al., 2001). Arterial compliance 
measurement was also conducted by monitoring internal pulsatile deformation in 

25 tissues surrounding the normal brachial artery (Mai and Insanna 2002). With some 
limits, these measurements have been correlated with clinical events including stroke 
(Buprez et al., 2001) and claudication symptoms (Taniwaki et al., 2001) in non- 
ESRD (End Stage Renal Disease) patients and adverse cardiovascular events in 
patients with ESRD (Blacher et al., 1998; Blacher et al., 1999), as well as length 

30 of time on dialysis (Luik et al., 1997). 
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One factor limiting the success of previously used methods is that 
arteries normally distended under physiologic pressure produce only small strain. 
The normal arterial wall, however, is a highly non-linear elastic medium, as 
illustrated by the solid curve in Figure 1 . The change of arterial elasticity due to 
5 intraluminal pressure was previously reported and analyzed over 40 years ago 
(Bergel 1961). Figure 1 qualitatively captures the essential feature of nonlinear 
arterial wall compliance. Under physiologic loading, the mean arterial pressure 
produces a high effective elastic modulus in the wall. Consequently, the arterial 
pressure pulse only creates small radial strain (Figure 1). 

10 Another factor limiting the success of previous methods is that 

properties of the vessel as a whole or in cross-section are measured. In the previous 
reports on the arterial compliance over a wide range of intraluminal pressure (Bank 
et al., 1999; Kaiser et al., 2001), the compliance was inferred from the geometrical 
changes such as artery diameter and lumen cross-section based on a numerical model 

15 (Langewouters* model; Langewouters et al., 1984). 

A phase-sensitive, two-dimensional speckle-tracking algorithm has 
been used by one of the co-inventors herein to determine displacements and strains 
(Lubinski et al., 1997). 

SUMMARY OF THE INVENTION 

20 An object of the present invention is to provide improved methods 

and systems for measuring mechanical property of a vascular wall and a method and 
system for determining health of a vascular structure wherein local deformation of 
a vessel wall resulting from physiologic pressures with altered transmural forces is 
measured. 

25 Using local measures of strain, highly localized measurements can 

be made. These local measurements can be used individually or compiled into a 
"map" or "image" of the mechanical properties of the vessel wall. Consequently, 
the methods and systems have the capability of being "high resolution." 
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In carrying out the above object and other objects of the present 
invention, a method for measuring a mechanical property of a vascular wall which 
deforms in response to a transmural force under usual physiologic pressures is 
provided. The method includes altering the transmural force to obtain an altered 
5 transmural force. The method also includes measuring local deformation of the 
vascular wall resulting from physiologic pressures with the altered transmural force, 
and determining a value for the mechanical property based on a measured amount 
of local deformation. 

The mechanical property may be a non-linear elastic property of the 

10 vascular wall. 

The step of measuring may include the step of non-invasively, 
ultrasonically imaging the vascular wall. 

The step of altering may include the step of reducing the transmural 
force to obtain a reduced transmural force. 

15 The step of reducing may include the step of applying an external 

pressure to the vascular wall. 

The external pressure may be substantially equal to a baseline internal 
pressure, and the vascular wall may deform by pulse pressure during a cardiac 
cycle. 

20 The step of reducing may include reducing an internal pressure to the 

vascular wall. 

The vascular wall may deform a relatively small amount in response 
to a transmural force under usual physiologic pressures and a relatively large 
amount in response to physiologic pressures with the altered transmural force. 
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The step of determining may include the step of directly estimating 
strain of the vascular wall. 

Further in carrying out the above object and other objects of the 
present invention, a method for measuring a mechanical property of a vascular wall 
5 is provided. The vascular wall is characterized by a relationship of arterial pressure 
versus strain that exhibits a relatively large slope under physiologic pressure caused 
by an arterial pressure pulse having a first mean arterial pressure and that exhibits 
a relatively small slope under physiologic pressure caused by an arterial pressure 
pulse having a second mean arterial pressure. The method includes altering the first 
10 mean arterial pressure to obtain the second mean arterial pressure. The method 
further includes measuring local deformation of the vascular wall at the second mean 
arterial pressure, and determining a value for the mechanical property based on the 
measured amount of local deformation. 

The step of measuring may include the step of non-invasively, 
15 ultrasonically imaging the vascular wall. 

The step of altering may include the step of reducing the first mean 
arterial pressure to obtain the second mean arterial pressure. 

The step of reducing may include the step of applying an external 
pressure to the vascular wall. 

20 The external pressure may be substantially equal to a baseline internal 

pressure, and the vascular wall may deform by pulse pressure during a cardiac 
cycle. 

The step of reducing may include reducing an internal pressure to the 

vascular wall. 

25 The step of determining may include the step of directly estimating 

strain of the vascular wall. 
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Still further in carrying out the above object and other objects of the 
present invention, a method is provided for determining health of a vascular 
structure which includes a vascular wall which deforms in response to a transmural 
force under usual physiologic pressures. The method includes altering the 
transmural force to obtain an altered transmural force. The method further includes 
measuring local deformation of the vascular wall resulting from physiologic 
pressures with the altered transmural force, and determining the health of the 
vascular structure based on the measured amount of local deformation. 

The step of measuring may include the step of ultrasonically imaging 
the vascular wall. 

The step of altering may include the step of reducing the transmural 
force to obtain a reduced transmural force. 

The step of reducing may include the step of applying an external 
pressure to the vascular wall. 

The external pressure may be substantially equal to a baseline internal 
pressure, and the vascular wall may deform by pulse pressure during a cardiac 
cycle. 

The step of reducing may include reducing an internal pressure to the 

vascular wall. 

The vascular wall may deform a relatively small amount in response 
to a transmural force under usual physiologic pressures and a relatively large 
amount in response to physiologic pressures with the altered transmural force. 

The step of determining may include the step of directly estimating 
strain of the vascular wall. 
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Yet still further in carrying out the above object and other objects of 
the present invention, a system for measuring a mechanical property of a vascular 
wall which deforms in response to a transmural force under usual physiologic 
pressures is provided. The system includes means for altering the transmural force 
5 to obtain an altered transmural force. The system further includes means for 
measuring local deformation of the vascular wall resulting from physiologic 
pressures with the altered transmural force, and means for determining a value for 
the mechanical property based on the measured amount of local deformation. 

The mechanical property may be a non-linear elastic property of the 

10 vascular wall. 

The means for measuring may include means for non-invasively, 
ultrasonically imaging the vascular wall. 

The means for altering may include means for reducing the 
transmural force to obtain a reduced transmural force. 

15 The means for reducing may include means for applying an external 

pressure to the vascular wall. 

The external pressure may be substantially equal to a baseline internal 
pressure, and the vascular wall may deform by pulse pressure during a cardiac 
cycle. 

20 The means for reducing may include means for reducing an internal 

pressure to the vascular wall. 

The vascular wall may deform a relatively small amount in response 
to a transmural force under usual physiologic pressures and a relatively large 
amount in response to physiologic pressures with the altered transmural force. 
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The means for determining may include means for directly estimating 
strain of the vascular wall. 

Still ftirther in carrying out the above object and other objects of the 
present invention, a system is provided for measuring a mechanical property of a 
5 vascular wall. The vascular wall is characterized by a relationship of arterial 
pressure versus strain that exhibits a relatively large slope under physiologic 
pressure caused by an arterial pressure pulse having a first mean arterial pressure 
and that exhibits a relatively small slope under physiologic pressure caused by an 
arterial pressure pulse having a second mean arterial pressure. The system includes 
10 means for altering the first mean arterial pressure to obtain the second mean arterial 
pressure. The system further includes means for measuring local deformation of the 
vascular wall at the second mean arterial pressure, and means for determining a 
value for the mechanical property based on the measured amount of local 
deformation. 

15 The means for measuring may include means for non-invasively, 

ultrasonically imaging the vascular wall. 

The means for altering may include means for reducing the first mean 
arterial pressure to obtain the second mean arterial pressure. 

The means for reducing may include means for applying an external 
20 pressure to the vascular wall. 

The external pressure may be substantially equal to a baseline internal 
pressure, and the vascular wall may deform by pulse pressure during a cardiac 
cycle. 

The means for reducing may include the means for reducing an 
25 internal pressure to the vascular wall. 
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The means for determining may include means for directly estimating 
strain of the vascular wall. 

Still further in carrying out the above object and other objects of the 
present invention, a system for determining health of a vascular structure including 
5 a vascular wall which deforms in response to a transmural force under usual 
physiologic pressures is provided. The system includes means for altering the 
transmural force to obtain an altered transmural force. The system further includes 
means for measuring local deformation of the vascular wall resulting from 
physiologic pressures with the altered transmural force, and means for determining 
10 the health of the vascular structure based on the measured amount of local 
deformation. 

The means for measuring may include means for non-invasively, 
ultrasonically imaging the vascular wall. 

The means of altering may include means for reducing the transmural 
15 force to obtain a reduced transmural force. 

The means for reducing may include means for applying an external 
pressure to the vascular wall. 

The external pressure may be substantially equal to a baseline internal 
pressure, and the vascular wall may deform by pulse pressure during a cardiac 
20 cycle. 

The means for reducing may include means for reducing an internal 
pressure to the vascular wall. 

The vascular wall may deform a relatively small amount in response 
to a transmural force under usual physiologic pressures and a relatively large 
25 amount in response to physiologic pressures with the altered transmural force. 
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The means for determining may include means for directly estimating 
strain of the vascular wall. 

The local deformation for the above methods and systems may be an 
intramural deformation. 

5 The above object and other objects, features, and advantages of the 

present invention are readily apparent from the following detailed description of the 
best mode for carrying out the invention when taken in connection with the 
accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

10 FIGURE 1 is a graph illustrating arterial elasticity characteristics at 

different preload; the solid curve represents normal tissue, and the dashed curve 
represents less compliant diseased tissue; 

FIGURE 2 is a schematic view of an experimental set-up for ex-vivo 
arterial elasticity measurements; 

15 FIGURE 3 illustrates representative image frames of a bovine carotid 

artery at different intraluminal pressures; from the top-left, pressure increases 
clockwise; the top three are in the low-preload region and the bottom three are in 
the high-preload region; each picture is separated by the same time interval; 
determining E from the arterial artery is very hard in the high-preload region 

20 because the arterial diameter hardly changes for a given pressure differential; 

FIGURE 4 is a graph of a pressure- strain curve for a bovine carotid 
artery; the lumen radius was 3.5 mm and the advential surface radius was 4.25 mm; 
the strain in the plot represents the radial normal strain in the middle layer of the 
artery wall as indicated by the dot in the inserted cross-sectional image; 
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FIGURE 5 are graphs which show accumulated radial displacement 
of two points on the artery wall of a healthy volunteer indicated by the two dots in 
the inserted cross-sectional image of a healthy volunteer; the dashed line represents 
the accumulated radial displacement of the point closer to the transducer and the 
5 solid line is for the point further away from the transducer; all lines are fit using 
cubic spline smoothing within an error of ± 0.01 mm. 

FIGURES 6a and 6b are graphs of radial strain and strain rate, 
respectively, of a healthy volunteer with continuous external compression; 

FIGURES 7a and 7b are graphs showing strain rate over one cardiac 
10 cycle; under physiologic pressure (i.e., Figure 7a), healthy (solid) and diseased 
(dashed) arteries are hard to distinguish; after pressure equalization (i.e., Figure 
7b), diseased (dashed) artery can be easily differentiated from healthy (solid) one; 

FIGURE 8a and 8b are schematic diagrams showing scanning of an 
upper arm without (i.e. , Figure 8a) and with (i.e., Figure 8b) a blood pressure cuff; 

15 while imaging the cross-section of the brachial artery and collecting ultrasound data 
frame-by-frame, surface compression was performed; the applied external pressure 
over the range of physiologic pressure is comparable to the pressure that would be 
generated in measuring a subject's blood pressure; both Figure 8a and 8b illustrate 
enough pressure equalization effect to bring the artery into the low preload region; 

20 external force can be monitored by a monometer attached to the cuff in Figure 8b; 

FIGURE 9 is a graph which illustrates radial normal strain referenced 
to the arterial geometry at maximum pulsation; the closed circles represent mean 
strain at each cardiac cycle, e ; and 

FIGURE 10 shows graphs which illustrate arterial elastic modulus 
25 reconstruction; elastic modulus represents *i[^J in semi-log scale and strain 
represents e ; open circles are represented moduli and solid line is the fit for the 
normal subject; open squares are reconstructed moduli and dashed line is the fit for 
the subject with known vascular disease. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

As previously mentioned, one of the significant obstacles of 
measuring mechanical properties of blood vessels is that vascular wall deformations 
under physiologic pressures only exhibit small deformations with the usual 

5 (physiologic) forces applied. By altering the forces across a vessel wall 
(transmural = across the vessel wall), the deformations induced from physiologic 
pressure across the arterial wall is altered with the methods and systems of the 
present invention. By measuring the local deformations, which may be intramural 
deformations, with altered transmural forces, new and detailed measurements can 

10 be made about the mechanical properties of the vessel. 

In other words,, one of the significant obstacles of measuring 
mechanical properties of blood vessels (arteries, veins, fistulae or other blood 
transporting structures hereafter referred to as "vessel" [noun] or "vascular" 
[adjective]) in the body (in-vivo) is that vascular wall deformations under 

15 physiologic pressures only exhibit small deformations with the usual (physiologic) 
forces (or pressures) applied. By altering the transmural forces across a vessel wall 
(transmural = across the vessel wall), the changes induced (deformation = changes 
in position of at least one location within the vessel wall) from physiologic pressure 
across the arterial wall is altered. By measuring the deformations with altered 

20 transmural forces, new and very detailed measurements can be made about the 
mechanical properties of the vessel. These measurements are derived from 
deformations induced by the altered forces. In particular, non-linear elastic 
properties of the vessel wall can be measured with great accuracy and precision if 
large wall deformations (e.g. , 20%) are induced. Non-linear elastic properties are 

25 directly related to many pathologies altering vascular compliance and, so, there 
should be numerous applications of a truly non-invasive method to monitor this 
physiologically significant parameter. 

The measurements include an parameter that can be derived from (1) 
altering the transmural forces (or pressures), and (2) measuring the local 
30 deformation that results from physiologic forces. These measurements include, but 
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are not limited to: displacement, strain, elastic moduli {e.g., Young's modulus), 
time derivatives {e.g., strain rate) or other mathematical operations or manipulations 
used to derive values from the measured local deformation. Consequently, an 
important aspect of the invention is measuring the local deformation of the vessel 
5 wall resulting from physiologic pressures with altered transmural force. 

The means of altering the transmural forces may be highly variable 
and include: external pressure applied to the vessel directly or applied to tissues 
overlying or surrounding the vessel such as, but not limited to: (1) manual 
compression {e.g., external compression with the hand, or the hand holding an 

10 object), (2) other methods of external compression such as a cuff inflated to apply 
pressure (such as, but not limited to a blood pressure cuff around an arm or leg), (3) 
internally "pressing" with a balloon or other internal structure that change geometry 
in such a way as to alter the transmural forces of a vessel, (4) maneuvers that alter 
forces such as straining against a closed airway (Valsalva maneuver) or blowing or 

15 sucking air against resistance, all of which change intra-thoracic pressure, alter 
blood flow in blood vessels with the thoracic cavity as well as blood entering and 
leaving the thoracic cavity, and thereby may change the transmural pressure (forces) 
of the vessel and allow the necessary measurements to be made, and (5) 
pharmacologic agents may be administered that alter blood pressure and thereby 

20 may be used to alter the transmural pressure, or (6) any other means of altering 
transmural vascular forces. 

The preferred means of measuring the local deformation uses data 
generated from ultrasound imaging. However, other means of measuring the 
deformation will allow calculation of the mechanical properties of the vessel. These 

25 other means of measuring the deformation may include, but are not limited to: (1) 
visual inspection, (2) manual measurement either directly or indirectly from images 
generated of the deformations from any imaging system, (3) sensing any location 
within the vessel wall with any method or device (such as, but not limited to, any 
wavelength electromagnetic radiation including, but not limited to, X-ray, CT, or 

30 measuring other physical parameters that allow the movement or position of the 
vessel wall to be measured such as MRI, or any imaging modality). All that is 

-14- 



UOM 0321 PUS 
2665 



important is to measure the location of the vessel wall to determine the local 
deformation from baseline (physiologic) pressures (forces) while the altered forces 
are generated. 

As described herein, the mean arterial pressure is lowered to reduce 
5 the preload so that the arterial pressure pulse creates much larger strain (Figure 1). 
The inherent elastic nonlinear ity of the arterial wall provides an opportunity to 
greatly expand the strain dynamic range by manipulating the mean arterial pressure. 

By lowering mean arterial pressure, it is much easier to differentiate 
diseased from normal arterial wall. Almost all arterial pathologies decrease 
10 compliance, as illustrated qualitatively by the dashed curve in Figure 1. At low 
arterial mean pressure (Figure 1), the difference in radial strain between normal and 
diseased arteries is much larger than at high mean pressure (Figure 1). An elasticity 
imaging procedure exploiting this relationship for highly sensitive characterization 
of arterial compliance is described herein. 

15 Arterial elasticity can be more accurately determined by measuring 

localized intramural strain. As described herein, the intramural radial normal strain 
is directly estimated using a phase-sensitive, two-dimensional speckle-tracking 
algorithm to determine displacements and strains (Lubinski et al., 1997). In a 
clinical setting, larger arterial strains with corresponding higher strain signal-to- 

20 noise ratio (SNR) are demonstrated using free-hand deformation to induce 
transmural pressure equalization and reduce preload. Strain and strain rate 
measurements at maximum pulsation correspond to the compliance of the artery 
under the same condition that blood pressure is taken with a blood pressure cuff. 
The feasibility of this technique is demonstrated using ex- and in-vivo measurements, 

25 and a straightforward elasticity reconstruction algorithm is presented to 
quantitatively assess the results. 

An example of preferred method of performing the measurement and 
resulting tenfold change in mechanical property (strain) measured is now described. 
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Peripheral Arterial Strain Imaging Using External Pressure Equalization 

Non-invasive peripheral arterial ultrasound strain imaging was 
performed while applying external pressure to induce changes in the baseline 
pressure difference across the arterial wall. By equalizing the baseline internal and 
5 external arterial pressure during the ultrasound measurement, increased arterial 
strains are induced by the relatively stable pulse pressure (pulse pressure = systolic 
pressure - diastolic pressure) during the cardiac cycle. The brachial and radial 
arteries of a 43 year old man were imaged with a 7.2 MHz linear ultrasound 
transducer while external deformations were continually increased over several 
10 cardiac cycles. External pressure was increased until the arteries collapsed during 
diastole, but distended during systole when the applied force exceeded the internal 
reference pressure of 80 mmHg (subject's diastolic pressure). Correlation-based, 
phase-sensitive, two-dimensional speckle-tracking algorithm was employed to 
calculate strain and strain rate from the internal displacement of the artery wall. 

15 Imaging of the arteries without external deformation resulted in 

measured strains up to 2% over the cardiac cycle. When applied pressure matched 
the internal baseline diastolic pressure of 80 mmHg, the strains increased by a factor 
of 10 with peak strains of 20% over the cardiac cycle. In addition, the peak strain 
rate under physiological conditions ranged from 0. 1 sec" 1 during diastole to -0.2 sec* 1 

20 during systole. After arterial pressure equalization, the peak strain rate increased 
to 1.0 sec 1 during diastole and -2.5 sec" 1 during systole, similar to the increase in 
peak strains. By applying external pressure, the pressure difference across the 
arterial wall at baseline (diastole) is reduced, while the pressure change from 
diastole to systole remains stable. As a result, the preload on the arterial wall can 

25 be decreased to near zero, leading to maximal strain during the cardiac cycle. By 
varying the external pressure, the range of measured strains vary over a cardiac 
cycle, and the non-linear properties of the arterial wall may be better characterized. 
Methods such as this that elicit the non-linear properties of the arterial wall could 
be used to better characterize vascular pathologies such as vessel hardening, 

30 neointimal hyperplasia, and vulnerable plaques. 
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As shown in Figure 1, if the mean arterial pressure can be lowered, 
then the arterial pressure pulse can create much larger strain. The inherent elastic 
non-linearity of the arterial wall provides an opportunity to greatly expand the strain 
dynamic range if the mean arterial pressure can be manipulated. 

5 By lowering mean arterial pressure, it is much easier to differentiate 

diseased from normal arterial wall. Almost all arterial pathologies decrease 
compliance, as illustrated by the dashed curve in Figure 1. At low arterial mean 
pressure, the difference in radial strain between normal and diseased artery is much 
larger than at high mean pressure. The focus here is to develop an elasticity 

10 imaging procedure exploiting this relationship for highly sensitive characterization 
of arterial compliance. 

Arterial elasticity can be more accurately determined by intramural 
strain. Preferably, the intramural radial normal strain is directly estimated. 

As described hereinbelow, in a clinical setting, larger arterial strains 
15 with corresponding higher strain signal-to-noise ratio (SNR) is demonstrated using 
free-hand deformation to induce transmural pressure equalization. Strain and strain 
rate measurements at maximum pulsation correspond to the compliance of the artery 
under the same condition that blood pressure is taken with a blood pressure cuff. 
The feasibility of this technique is demonstrated using ex- and in-vivo measurements. 

20 Methods 



Like any other tissue, arteries exhibit non-linear elasticity (Fung 
1993). To demonstrate and quantify arterial non-linearity with respect to internal 
loading, a controlled experimental protocol was designed. The intramural strain of 
an ex-vivo bovine artery was measured when the internal pressure was increased by 
25 a fluid-filled syringe pump. 

Experimental Set- Up . A closed-loop compression system was 
designed to pressurize an artery sample while simultaneously scanning with 
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ultrasound. A programmable commercial syringe pump (Cole-Parmer) serves as a 
pressure source. An acoustic window is designed to hold ail arterial sample between 
inlet and outlet ports. The outlet is sealed so that' the internal pressure develops 
while the syringe pump compresses. A pressure gauge is placed between the 

5 syringe and artery sample close to the artery to measure intraluminal pressure. The 
acoustic window is placed in a water tank with anechoic material at the bottom to 
suppress possible reverberation, and the tank is placed underneath an ultrasound 
transducer positioning device. A PC -based RF data acquisition system is connected 
to a commercial ultrasound scanner (Siemens Elegra). A block diagram is presented 

10 in Figure 2. 

Bovine Artery ex-vivo . A 50 mm-long bovine carotid artery segment 
preserved in 30% ethanol (Artegraft, Brunswick, NJ) was placed in the middle of 
the acoustic window connected to the flow path filled with degassed water. A 
commercial (Cole-Parmer) syringe pump was programmed to pump water at a fixed 

15 rate over a fixed period (70.5 ml/min. for 13 seconds) to build intraluminal pressure 
to 120 mmHg. While the artery distended from the resting position, a 12.0 MHz 
linear ultrasound array connected to a commercial ultrasound scanner (Siemens 
Elegra) imaged the arterial cross-section at a rate of 22 frames per second for 13 
seconds. RF data from every frame in the sequence were captured. The 

20 intraluminal pressure over time was also recorded. Data were subsequently 
processed using a phase-sensitive, two-dimensional speckle-tracking algorithm to 
determine displacements and strains (Lubinski et aL, 1997). Correlation-based 
algorithms were used to track internal displacements. Frame-to-frame lateral and 
axial displacements were estimated from the position of the maximum correlation 

25 coefficient, where the correlation kernel size equaled the speckle spot for optimal 
strain estimation and axial displacements were refined using the phase zero-crossing 
of the complex correlation function. Frame-to-frame displacement estimates were 
integrated from and registered to the initial coordinate system (i.e., Lagrangian 
presentation). Spatial derivatives of the displacements were computed in one region 

30 of the artery to estimate the radial normal strain (i.e., the radial derivative of the 
radial displacement). As described hereinbelow, the radial normal strain is called 
simply the strain, where appropriate. 
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Human Artery in-vivo . Two subjects were tested. The first was a 43 
year old healthy male volunteer and the second was a 48 year old male with ESRD 
secondary to diabetes mellitus, on hemodialysis, and a history of peripheral vascular 
occlusive disease, with prior right below the knee amputation. 

5 A 7.2 MHz linear array was used with continuous freehand 

compression performed on the surface of the right upper arm close to the brachial 
artery. While imaging the cross-section of the brachial artery at a rate of 107 
frames per second and collecting ultrasound data frame-by-frame, surface 
compression was performed by the investigators. The applied external force 

10 produces internal pressure comparable to that generated in measuring a subject's 
blood pressure. The compression was increased until brachial artery pressure 
exceeded diastolic pressure, as evidenced by viewing B-scan images. Collected RF 
data were processed off-line in the same way as described above. 

Results 

15 Bovine Artery ex-vivo . The accumulated displacement and axial 

normal strain within the artery wall were referenced to a frame where the artery was 
at rest. Representative image frames over the full pressure range are presented in 
Figure 3. From the resting position of the bovine artery at 0 mmHg internal 
pressure, the syringe pump pressurized the lumen to 120 mmHg. The pressure is 

20 plotted in Figure 4 versus axial normal strain at the arterial bottom wall, 
approximating the radial strain. The strain initially increases rapidly with pressure 
and then hardly changes. Assuming incompressibility and plane strain in the cross- 
section of the artery, the radial normal strain, e n on the artery wall can be expressed 
as a function of internal pressure, p h outside pressure, p 0 and Young's modulus, E 

25 (Appendix A and Figure 1 1): 

e -e, - ■ \ ah \ 2 f a<r<b, (1) 
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where a is the lumen radius, b is outer radius of the artery, and r is the strain 
measurement point. Using this relation, E can be estimated as a function of e. In 
this experiment, outside pressure, p 0 can be ignored because the artery sits less than 
2.5 cm from the surface of the water. As illustrated in Figure 4, the elastic modulus 
5 at low preload (E=39kPa at 40% preload, representing an intramural pressure of 
about 3kPa) is over an order of magnitude smaller than that at physiologic preload 
(E=454kPa at 60% preload, representing an intramural pressure of about lOkPa). 
This implies that the artery in the physiologic region produces only 1/10 of the 
strain in the low-preload (i.e., after pressure equalization) region for the same 

10 pressure differential. Because it is difficult to distinguish diseased from normal 
arteries with small strain under physiologic pressure, it is important to bring the 
artery into the low-preload region to better determine its elasticity through strain 
measurement. By applying external force opposite to that generated by the mean 
arterial pressure, higher strains may enable easier differentiation of diseased arteries 

15 from normal ones. 

Human Artery in -vivo . The accumulated radial displacement of the 
brachial artery wall was estimated relative to the original frame, as illustrated in 
Figure 5, where results are compared from two points in the wall of the normal 
subject separated by 0.2 mm. The overall slope of the average displacement signals 
20 compression by the transducer. Cyclic displacement represents deformation by the 
pulse pressure. The point on the bottom wall closer to the transducer displaced 
more than the point further away from the transducer at systole and both points 
come back to the resting position at diastole. This shows that the artery wall 
compresses maximally at systole and returns to the resting position at diastole. 

25 Based on the displacement information depicted in Figure 5, radial 

normal strain and strain rate were estimated. When the applied pressure matched 
the internal baseline diastolic pressure of 80 mmHg, strains increased by a factor of 
10 with peak strains of 20% over a cardiac cycle (Figure 6a). In addition, the strain 
rate computed from the accumulated displacements increased over several cardiac 

30 cycles. The peak strain rate under physiologic conditions ranged from 0.1 Hz 
during diastole to -0.2 Hz during systole. After arterial pressure equalization, the 
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peak strain rate increased to 1.0 Hz during diastole and -2.5 Hz during systole, 
similar to the increase in peak strains (Figure 6b). By equalizing the baseline 
arterial pressure to approximate the diastolic pressure, the preload on the arterial 
wall decreases to zero, resulting in maximal strain. The elastic properties of the 
arterial wall can be better characterized with intramural strain measurements 
extending over a large preload range. 

The radial normal strain and strain rate of the diseased subject are 
compared with those of the normal, over approximately one cardiac cycle in Figures 
7a and 7b. Figure 7a illustrates differences at high preload (i.e., physiologic 
condition), and Figure 7b at an intraluminal pressure nearly equalized by the applied 
force (i.e., low preload). Overall, the artery of the diseased subject produces 
smaller strain and strain rate than the artery of the normal. Pulse pressure of each 
subject was recorded by measuring blood pressure before and after the ultrasound 
scan. The blood pressure of the diseased subject was 160/85 mmHg while the 
normal subject measured 114/78 mmHg. Noting that the pulse pressure (Ap = 75 
mmHg) of the diseased subject is twice as big as the normal subject (Ap = 36 
mmHg), the strain and strain rate contrast will be doubled after normalization by the 
pulse pressure. These results demonstrate that differences in elastic properties 
between the two subjects become more pronounced when measured in the low- 
preload region by this pressure equalization procedure. 

Elastic Modulus Reconstruction . For an isolated artery, Equation (1) 
can be used to reconstruct the arterial elastic modulus over the entire strain range, 
as illustrated in Figure 3. However, in-vivo measurements cannot use this simple 
formulation since the artery is attached to surrounding tissue. If surrounding tissue 
can be modeled as a continuous medium with elastic modulus E 2 , and the artery wall 
is considered homogeneous with an elastic modulus Ej, then the radial strain can be 
written as (Appendix B and Figure 12): 



e = e = 



-la 2 b 2 Pi 
2(b 2 -a l )r l 



(2) 
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To reconstruct the modulus, this equation must be inverted, 



-3a Vp, 
2(b 2 -a 2 )r 2 



e 



(3) 



As a first step in solving this equation, a, b, and r must be estimated 
from B-scan images. In this study, the constants were computed by hand, but in . 
5 real-time clinical operation it is very feasible to design automatic lumen detecting 
algorithms to define both intimal and advential boundaries (i.e., a and b). To 
determine a, b, and r, the averaging procedure was used in this study because 
artery was not perfectly circular over the entire pressure equalization procedure. 
This averaging should not introduce significant error into the estimated modulus as 
10 long as a and b change at the same rate for the reasonably small deformation. 

Given a and b and the coordinates of the strain measurement position, 
the radius r can also be captured automatically. Consequently, Equation (3) can be 
written as: 

E X +K 2 E 2 =*i[|j]> ( 4 > 

_ 3 2.2 2 

15 where K x = — - = — -, K 2 =— — are geometric factors computed from B-scan 

2(b —a )r (b —a ) 

images. Ap is pulse pressure and Ae is inter-cardiac strain (i.e., change in strain 
from systole to diastole). The strain (Figure 6a) for the normal subject referenced 
to the arterial geometry at high preload (under physiologic pressure) is converted 
into the strain referenced to the arterial geometry at low preload (after pressure 
20 equalization) to conform to the normal geometry used to compute non-linear elastic 
parameters. That is, the reference frame must be converted to the one closest to the 
undistended arterial geometry to present Lagrangian strain for nonlinear analysis 
similar to that presented in Figure 4 for the isolated artery. This means frame-to- 
frame displacements must be accumulated "backwards" from the last frame to the 
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first, yielding the final simple expression relating the strain relative to the first frame 
( e° for the frame) to the strain relative to the last frame (e, for the i th frame): 



£ _ £ N-i+\ £ N 



(5) 



where e° N is the maximum original strain at the pressure equalization frame (Figure 
5 9). 

Inter-cardiac strain of varying amplitude developed over each cardiac 
cycle Aeand corresponding mean strain e was calculated from Figure 9. With 
intimal and advential boundaries, a and b, and the coordinates of the strain 
measurement position, r estimated from B-scan images, ^j^-j from Equation (4) 

10 is plotted on a semi-log scale versus mean strain at each cardiac cycle, e (Figure 
10). The pulse pressure measured before and after the experiment, Ap, was 
assumed to be constant over cardiac cycles. The reconstructed elastic moduli were 
fit to a straight line. If the elastic modulus of the surrounding muscle, E 2 , can be 
considered small compared to the arterial elastic modulus, the intercept will 

15 determine the undistended (i.e., zero preload) in-vivo arterial elastic modulus. 
Otherwise, the elastic modulus of surrounding muscle must be measured to correctly 
reconstruct the arterial elastic modulus, which is overestimated by KJE 2 . However, 
if change in elastic modulus with preload is small for surrounding tissue compared 
to that of the artery wall, then the slope of this curve (i.e. , nonlinear coefficient for 

20 a purely exponential model) will be correctly reconstructed independent of the 
elastic modulus of the surrounding medium. For the normal subject, the intercept 
ranges from 14.7 kPa to 16.5 kPa, and the slope is 2.9 within an error of ± 0.1. 
The intercept ranges from 153.2 kPa to 193.7 kPa, and the slope is 4.0 within an 
error of ± 0.6 for the subject with known vascular disease. The undistended elastic 

25 modulus of the diseased subject is over ten times that of the normal subject. The 
ratio of the elastic modulus of the diseased artery over the normal artery will only 
increase when accounting for the overestimation by K^ 2 from the surrounding 
tissue, which is assumed to be small compared to £*,. Since the surrounding tissue 
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modulus may be comparable between subjects, the overestimation for the diseased 
subject is much less pronounced than that for the normal subject. In addition, the 
nonlinear coefficient (slope in Figure 10) is more increased for the diseased subject. 

The nonlinear coefficient can serve as a strong indicator of arterial 
5 stiffness. To estimate this parameter with optimal accuracy, all image frames from 
high preload to low must be used. One limitation of the present study was that all 
image frames for the diseased subject were not of high enough quality to contribute 
to this accumulation. Consequently, high precision intramural strain was computed 
over several non-consecutive cardiac cycles in the image sequence. An absolute 
10 geometric reference was established between isolated segments by tracking changes 
in arterial wall thickness over the entire sequence using B-scan images. These lower 
precision measurements only provided the geometrical reference for high precision 
intramural strain measurements. Nevertheless, the results presented in Figure 10 
show that good results can be obtained from intramural measurement even when all 
15 image frames are not of sufficient quality for high precision phase-sensitive speckle 
tracking. 

Discussion and Summary 

The intramural strain range in peripheral arteries produced by the 
pulse pressure can be significantly extended by simply applying pressure comparable 

20 to a subject's blood pressure. Intramural strain can be monitored directly with high 
precision using phase-sensitive ultrasonic speckle tracking algorithms developed for 
elasticity imaging. By combining pressure equalization with phase-sensitive speckle 
tracking, new diagnostic information may be gathered about the non-linear elastic 
properties of the arterial wall. As demonstrated experimentally herein, the radial 

25 strain and strain rate increased ten-fold in a healthy artery when mean arterial 
pressure is reduced from physiologic levels. The deformation in a diseased artery, 
however, changed comparatively little as the pressure was equalized. Consequently, 
a diseased artery was easily differentiated from normal simply by observing radial 
normal strain and strain rate during the compressed phase of the examination. 

30 These very preliminary ex- and in-vivo results suggest that even small changes in 
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arterial stiffness accompanying vascular disease may be sensitively monitored with 
elasticity imaging. 

In addition to qualitative assessment of vascular compliance, the non- 
linear elastic modulus of the vascular wall can be quantitatively estimated using a 
5 simple reconstruction procedure. If surrounding tissue can be modeled as a 
continuous medium with elastic modulus E 2 , the elastic modulus Ej of the artery 
wall can be reconstructed using Equation (4). Within an offset proportional to E 2 , 
it is possible to reconstruct the arterial elastic modulus as a function of mean arterial 
strain from the ratio j^-j at the following different levels of sophistication: 

10 > assume a fixed geometric scale factor K l for all subjects; 

► automatically determine a and r from B-scan images using lumen 
boundary detection algorithms assuming a fixed wall thickness for all 
subjects to compute K t \ 

► automatically determine a t b, and r from B-scan images to compute 
15 Kj. 

This procedure overestimates the arterial elastic modulus by as much as This 
should not pose a practical problem at high preloads, but it could be a source of 
error at low preloads (i.e. , pressure equalized) where K^ 2 may not be much smaller 
than E } . The offset may only be a practical issue for low arterial moduli where it 
20 can limit the detection of subtle changes in arterial compliance such as in normal or 
near normal arteries. 

Both ex- and in-vivo measurements presented herein, as well as a 
large body of previous literature, suggest that the non-linear change in arterial 
elastic modulus with preload can be modeled as an exponential function. 
25 Consequently, a simple linear least squares fit to the natural log of the estimated 
elastic modulus as a function of preload can fully characterize the vessel wall's non- 
linear mechanical properties. A major advantage of this fit procedure is that only 
a few points are needed over a limited range of preloads to estimate the elastic 
modulus of the undistended artery. This may be very valuable in applications such 

-25- 



UOM 0321 PUS 
2665 



as assessment of carotid compliance where it may be difficult to equalize the 
pressure all the way to the diastolic limit. Again, the elastic properties of the 
background medium will influence the fit, but they should not significantly alter "the 
results except in the small preload limit of highly compliant arteries. In any event, 
5 both the intercept (i.e. , the elastic modulus of the undistended artery) and the non- 
linear parameter (i.e. , the slope of the fit) can assess the vascular compliance. 
Because of the large non-linear parameter in arterial tissue compared to most soft 
tissues, the slope should not be greatly affected by the surrounding medium. This 
hypothesis, as well as the predictive value of each parameter, can be tested in 
10 controlled ex-vivo studies. 

External force measurement provides additional information about the 
elastic modulus of the surrounding muscle to further refine the reconstruction 
procedure. A commercial blood pressure cuff has been modified to have an acoustic 
window through which an ultrasound scan can be performed. A monometer 
15 attached to the cuff monitors the external force (Figures 8a and 8b). Force is 
recorded and input to a modified reconstruction algorithm accounting for the finite 
elastic modulus of surrounding tissue. This procedure can also be tested in 
controlled ex-vivo studies. 

Reconstruction procedures presented herein focus on the static elastic 
20 properties of the arterial wall. It is well known that the arterial wall is a viscoelastic 
medium (Hardung, 1962). Consequently, additional information can be obtained by 
comparing strain rate measurements with the arterial pressure pulse to derive time 
constants related to viscoelastic parameters. 

Assessing arterial elasticity has many important clinical applications. 

25 This method allows localized assessment of vascular elasticity that may reflect the 
degree of both local and general vascular disease. It may be useful in pre-operative 
assessment for certain vascular surgery procedures, since the elastic properties of 
the vessel may reflect the capacity of the artery to remodel, influencing clinical 
outcomes. For example, in surgically creating an arterial- venous anastomosis in 

30 hemodialysis fistula creation, the artery dilates to create a manifold increase in 
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volume flow through the fistula to accommodate hemodialysis (Konner et aL, 2003). 
Inelastic, diseased arteries, so prevalent in end-stage renal disease, may greatly 
influence the outcome of the procedure (Konner et al., 2003). Assessing the 
elasticity of arteries preoperatively may favorably influence site selection, prevent 
the development of peripheral ischemia and improve clinical outcomes. The ease 
of collecting data reliably, such as with a modified blood pressure cuff (Figure 8b), 
will be important in assessing the utility of this method in the clinical setting. 
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APPENDIX A 



Referring now to Figure 11, consider an elastic circular cylinder 



subject to hydrostatic pressure. The cylinder is long enough compared to its cross- 
sectional area to assume a plane strain state. Assuming isotrophy in the axial 
5 direction and homogeneity in the radial direction, only normal stress exists. 
Considering axial symmetry, radial displacement depends only on the distance r 
from the center of the cylinder: 



where A and B are constants related to the material characteristics and geometry. 
10 The corresponding radial and tangential strains are: 



u r {r) = + flr, 



(A-l) 



r 




(A-2) 




-4 



(A-3) 



. Stress-strain relations for this problem are: 




(A-4) 



15 




(A-5) 




(A-6) 



where o, is the /'* component of the stress tensor and v is Poisson's ratio. 
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Since e z = 0 for a plane strain case, Equation (A-6) becomes: 



a z =v(cr r +a e ): (A-7) 
With this expression, Equations (A-5) and (A-6) can be rewritten as: 



e r (r) =-^(1 +v)[(l -v)a r -vcx e ] (A-8) 



e e (r) = -Jrd + v)[(l -v)<r 0 -v(T f ], (A-9) 



Combining Equations (A-8) and (A-9) to eliminate a e yields: 



(1 -v)e, + =-Jr(l + v)(l -2v)a> . (A- 10) 

E 



Substituting Equations (A-2) and (A-3) and applying boundary conditions at inner 
(a r = - Poj a t r = where p Q is the external pressure) surfaces lead to: 



(\-2v)Ar + B = — ^(i+v)(l-2v) (A-ll) 
a" E 



(1 -2v) 4" + * = — +v)(l -2v). (A-12) 



Combining Equations (A-ll) and (A-12) determines the unknown coefficients, A and 
B. The radial strain of interest herein can be expressed as: 
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_ { Po -p,) 0 V(l-v) {b 2 p 0 -a 2 P| .)(H-v)(l-2v) 



E{b 2 -a 2 )r 2 E(b 2 -a 2 ) 



Assuming incompressibility, v = 0.5, Equation (A-13) reduces to: 



- p ' )a 2 t a<r<b , (A . 14) 
2E(b 2 -a 2 )r 2 
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APPENDIX B 

Referring now to Figure 12, consider the same cylinder as in 
Appendix A surrounded by a homogeneous material which has different elastic 
modulus. Employing the same assumptions as in Appendix A, displacements and 
5 resulting strains in media I and II can be expressed as follows: 



«/.r(0=- — + Br (B-l) 
r 



or r 1 



£,,e( r ) = ! y= ~ (B-3) 



II: u llr (r) =-— + Dr (B-4) 

r 



10 €nA r)=^ = ^- + D (B-5) 

or r 



e lt e(r) = y= -jf+D. (B-6) 



Note that D=0 to satisfy the boundary condition at infinity, i.e. , no displacement. 
Using the same procedure as in Appendix A, radial stresses in medium I and II can 
be expressed as follows: 



15 I: ^ (r) = (l + v)0-2v) t(1 - v)£ ^- + V£ ^ ] (B " 7) 
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H: a "-' ir) = ( l + v)U-2v) [(1 " V)E "r + v£ »^ 



where £ y is the arterial elastic modulus and E 2 is the modulus of the surrounding 
medium. Applying a boundary condition at the inner surface (a / r = - p if at r = a), 
and two boundary conditions at the outer surface (a I/r = - p 0 , a IIr = o ln at r = b) 
5 yields: 

(l-2v)A + B = ^£L ( i +v)( i_2v) (B-9) 



-— + tB = _£ (B-10) 
ft b 



(l-2v)£ 2c= (l-2v)£ lA + 
f> 2 i 2 



Combining Equations (B-9), (B-10), and (B-ll) determines A, B, and C: 



A p,(l+v)fl 2 fc 2 [E, +(l-2v)E 2 ] (B12) 

£ 1 {(£,-£ 2 )a 2 -[£,+(l-2v)£ 2 ]t 2 } 



10 



B -p,.(l+v)(l-2v) fl 2 (£,-£ 2 ) (B13) 

£,{(£, -£ 2 )a 2 -[£, +(l-2v)£ 2 ]t 2 } 



C 2p,-(l+v)(l-v) a V£, (B _ 14) 

£,{(£, -£,) fl 2 -(£,+(l-2v)£ 2 ]ft 2 )' 
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Assuming incompressibility, v = 0.5, radial strain in medium I can be reduced to: 



-3aVp,- 
2(fc 2 -a 2 )r 2 



(B-15) 



While embodiments of the invention have been illustrated and 
described, it is not intended that these embodiments illustrate and describe all 
possible forms of the invention. Rather, the words used in the specification are 
words of description rather than limitation, and it is understood that various changes 
may be made without departing from the spirit and scope of the invention. 
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